Energy pumping in a quantum nanoelectromechanical system 
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Fully quantized mechanical motion of a single-level quantum dot coupled to two voltage biased 
electronic leads is studied. It is found that there are two different regimes depending on the applied 
voltage. If the bias voltage is below a certain threshold (which depends on the energy of the 
vibrational quanta) the mechanical subsystem is characterized by a low level of excitation. Above a 
threshold the energy accumulated in the mechanical degree of freedom dramatically increases. The 
distribution function for the energy level population and the current through the system in this 
regime is calculated. 



During the past few years experimental methods of 
physics has seen an advancing capability to manufac- 
ture smaller and smaller structures and devices. This has 
lead to many new interesting investigations of nanoscale 
physics. Examples include, for instance.observation of 
the Kondo effect in single-atom j unctions [ij], manufactur- 
ing of single- molecular transistors 2], and so on. There 
has also been a great interest in the promising field of 
molecular electronics 3J. One of the main features of 
the conducting nanoscale composite systems is its sus- 
ceptibility to significant mechanical deformations. This 
results from the fact that on the nanoscale level the me- 
chanical forces controlling the structure of the system are 
of the same order of magnitude as the capacitive electro- 
static forces governed by charge distributions. This cir- 
cumstance is of the utmost importance in the so called 
electromechanical single-electron transistor (EM-SET), 
which has been in focus of recent research. The EM-SET 
is basically a double junction system where the additional 
(mechanical) degree of freedom, describing the relative 
position of the central island, significantly influences the 
electronic transport. Experimental work in relation to 
EM-SET structures range from the macroscopic |5| to the 
micrometer scale H H 13 and down to the nanometer 
scaled . Various aspects of electronic transport in such 
systems have been theoretically investigated in a series 
of articles [TJ El El El El El El El 111 Ei 

In Ref. and Ref. it was, among other things, 
shown that coupling the mechanical degree of freedom 
of an EM-SET to the nonequilibrium bath of electrons 
constituted by the biased leads, can lead to dynamical 
self excitations of the mechanical subsystem and as a re- 
sult bring the EM-SET to the shuttle regime of charge 
transfer. This phenomena is usually referred to as a shut- 
tle instability. In these papers the grain dynamics are 
treated classically and the key issue is that the charge of 
the grain, q(t), is correlated with it s velocity, x(t), in a 
way so that the time average, q{t)x ^ 0. 

Decreasing the size of the central island in the EM- 
SET structure to the nanoscale level results in the quan- 
tization of its mechanical motion. Charge transfer in 
this regime was studied theoretically in Ref. El How- 



ever, the strong additional dissipation in the mechanical 
subsystem suggested in this paper keeps the mechanical 
subsystem in the vicinity of its ground state and prevents 
the developing of the mechanical instability. The aim of 
our paper is to investigate the behavior of the EM-SET 
system in the quantum regime when its interaction with 
the external thermodynamic environment generating ad- 
ditional dissipation processes can be partly ignored in 
such a manner that the mechanical instability becomes 
possible. We will show that in this case at relatively 
low bias voltages, intrinsic dissipation processes bring the 
mechanical subsystem to the vicinity of the ground state. 
But if the bias voltage exceeds some threshold value, the 
energy of the mechanical subsystem, initially located in 
the vicinity of the ground state, starts to increase expo- 
nentially. We have found that intrinsic processes alone 
saturate this energy growth at some level of excitation. 
The distribution function for the energy level population 
and the current through the system in this regime is cal- 
culated. 

We will consider a model EM-SET structure consist- 
ing of a one level quantum dot situated between two leads 
(see Fig. ^| . To describe such a system we use the Hamil- 
tonian 

H = ^%«L«M + (£o-^)c , C +^ 2 

A;, a 

+ ima)„X 2 + T k:a {X)[a\ a c + c 1 a k , a }. (1) 

The first term describes electronic states with energies 
Ek, a and where a' ka (a^ iQ ) are creation (annihilation) 
operators for these noninteracting electrons with momen- 
tum k in the left (a = L) or right (a — R) lead. The 
second term describes the interaction of the electronic 
level on the dot with the electric field so that (c) 
is the creation (annihilation) operator for the dot level 
electrons and Eq is the dot energy level. The scalar 
D represents the strength of the Coulomb force acting 
on a charged grain, X is the position operator, and 
xq = mwo is the harmonic oscillator length scale 
for an oscillator with mass m and angular frequency u>o- 
The third and fourth terms describe the center of mass 



movement of the dot in a harmonic oscillator potential ( a ) 
so that P is the center of mass momentum operator. 
The last term is the tunneling interaction between the 
lead states and the dot level and T ka (X) is the tunnel- 
ing coupling strength. We will consider the case when 
the tunneling coupling depends exponentially on the po- 
sition operator X, i.e. Tk t n(X) = T.Rexp{X/A} and 
T k>L (X) = T L exp{-X/A}, where T R and T L are con- 
stants and A is the tunneling length. 

To introduce a connection to the quantified vibrational 
states of the oscillator we perform a unitary transforma- 
tion of the Hamiltonian Q so that UHW = H, where 3 ) 
U = e^ pd ° c c . In this paper we consider the situation 
when H has the most symmetric form: 

H = E k:a a\ a a k ^ a + E Jc + huj (tfb + ^ j 
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Here w (b) is a bosonic creation (annihilation) opera- 
tor for the vibronic degree of freedom, and the dimen- 
sionless parameters x± = 1/\/2{xq/ K ± do/x$) (where 
do = D / (xoitiloq)) characterize the strength of the elec- 
tromechanical coupling. Furthermore, To is the renor- 
malized tunneling coupling constant, and Eq — Eo — 
Ddo/(2xo + rauj^d\l2) is the shifted dot level. For sim- 
plicity, but without loss of generality, we choose Eq equal 
to the chemical potential of the leads at zero bias voltage. 

First let us study the situation when the mechanical 
subsystem is characterized by a low level of excitation. 
We will consider the case of small electromechanical cou- 
pling. This means that the dimcnsionless parameters 
x± <C 1 and that only elastic electronic transitions and 
transitions accompanied by emission or absorption of a 
single vibronic quantum (single-vibronic processes) are 
important. If the applied voltage is smaller than 2huJo/e 
and the temperature is equal to zero, the six allowed 
transitions of this type are the ones described in Fig. HJi. 
Here we see that only elastic tunneling processes and tun- 
neling processes in which the vibronic degree of freedom 
absorbs one vibronic quantum are allowed and as a re- 
sult the rate equation for the distribution function of the 
energy level population P(n, t) has the form: 

T- l d t P{n,t) = P(n,t) + (x 2 + + x 2 _)(n + l)P(n+l,t) 
- (1 + n{x\ + x 2 _)) P{n,t), 

where T = 2ttTq v/h and v is the density of states in the 
leads. 

It is straightforward to solve these equations and find 
that the solution exponentially fast approaches the stable 
solution P(0) = 1 and P{n) = for all n > 0. As a 
result, the dimensionless average extra energy excited in 



FIG. 1: (a) Model system consisting of a one level quantum 
dot placed between two leads. The level of the dot equals 
the chemical potential of the leads and a bias voltage of V 
is applied between the leads. The center of mass movement 
of the dot is in a harmonic oscillator potential with the vi- 
brational quanta Hujo- The applied bias voltage is such that 
eV/2 < huJo- (b) Same as (a) but with the applied bias voltage 
larger than 2huo/e. 



the vibronic subsystem, 



(3) 



n=0 



goes to 0. 

If the applied bias voltage is increased above the 
threshold value V c = 27j,wo/e we instead get the allowed 
transitions described in Fig. i.e. two absorption pro- 
cesses has changed into emission processes where the en- 
ergy quantum Hloq is transferred to the vibronic degree of 
freedom. These transitions lead to the following equation 
for P(n): 

Y- l d t P{n,t) = -(l + x 2 _ + {n+l)x 2 + )P(n,t) 

+ l)P(n + !,*) + P(n, t) + x\nP{n - 



One can find from this equation that the time evolution 
of the exited energy is given by the formula: 



E(t) = 



i 



(5) 
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FIG. 2: Illustration of the second order case where elastic 
tunneling and inelastic tunneling exchanging two or less vi- 
brational quanta are included. The level of the dot is equal 
to the chemical potential and the bias voltage is set so that 
2huj /e < V < 4?kj /e. 



i.e. energy is continuously pumped into the mechanical 
subsystem, which is strong evidence that the low exited 
regime is unstable if the bias voltage exceeds the critical 
value V c . Furthermore, it is necessary to remark here 
that for this case we thus have a linear increase in the 
energy as a function of time even when x 2 , approaches 
x 2 _. 

As the excitation of the vibronic subsystem increases 
multi-vibronic processes become important. They give 
rise to an additional dissipation which saturates the en- 
ergy growth induced by the single- vibronic processes. As 
a result the system comes to a stationary regime which 
is characterized by a significant level of excitation of the 
vibronic subsystem. To demonstrate this we will now 
expand our analysis by taking into account electronic 
transitions accompanied by the emission or absorption 
of two vibronic quanta (two- vibronic processes). To de- 
scribe such transitions one has to take into account sec- 
ond order terms in w and b in the tunneling part of the 
Hamiltonian @. As illustrated in Fig. |2]these terms will 
generate four processes in which two vibrational quanta 
are absorbed by the electron during the tunneling event. 
There is also a renormalization of the elastic channel com- 
ing from the inclusion of these terms. Now the equation 
for the distribution function of the energy level popula- 
tion has the form: 



Y- l d t P{n,t) = nP{n- l,t) 

- en 2 + (a - e + l)n + 1 
+a(n + l)P(n + l,t) 



P(n,t) 



+e(n + l)(n + 2)P(n + 2,t) = 0, (6) 
where we have introduced the constants e = (4 + 



xi)/(4a;^) and a = x 2 _/x 2 + . 

To find the stationary solution of this equation we in- 
troduce the generating function: 

oo 

P(z) = ^>"P(n), 
n=Q 

where z is a complex number inside the unit circle. 
Rewriting Eq. Hjjwe find the equation for V(z) 

e(z + l)d 2 V{z) + (a- z)d z V(z) - V(z) = 0. 

The solution to this equation is 

v{z) = e ~fi dz '^4 

\J Z0 e(*'+l)J 

where C\ and zq are constants. Since the proba- 
bilities P(n) are positive and normalized, the sum 
^^ (— l)™P(n) = V(z — —1) converges absolutely. 
This is true only for zq = — 1. The second constant 
C\ can be determined from the normalization condition 
V{z = 1) = 1 to be Ci = e2T/ dxe^- x ^ e (x + l) 7 " 1 , 
where we have introduced the constant 7 = (a — e+l)/e. 
Therefore the final expression for V{z) is 



V(z) 



Ci e" 



e (z + 1)t y_ 



dz'e-^{z' + 1)T- 1 . (7) 



We can now calculate the average energy excited in 
the harmonic oscillator, which is just d z V{z) calculated 
at z = 1, 



E=^{2 + C 1 -e 1 ). 



(8) 



To see how the energy pumped into the harmonic oscil- 
lator affects the charge transport we calculate the current 
/ through the system in units of eT. For voltages below 
V c the current is only mediated by the elastic channel 
and is thus / — 1/2. 

For voltages in the range 2hujo/e < V < Ahuio/e the 
current can be calculated to leading order in e as 

I = - 4) + x\x 2 _^j V"(z = 1) + (x + - x_f V'(z = 1) + l.( 

In Fig. we have chosen a set of numerical values and 
plotted (solid line) the calculated current as a function 
of the bias voltage. For comparison we have also plotted 
(dashed line) the current as given in the high dissipation 
limit where the harmonic oscillator goes to the ground 
state between tunneling events. It is clear that the cur- 
rent in the regime characterized by a high level of exci- 
tation is much greater than the one in the regime of low 
level of excitation. 
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FIG. 3: The current as plotted as a function of the bias voltage 
(solid line). The current makes a jump as the non-elastic 
channel is opened at V = 2Tiu>o/e. We have also plotted 
(dashed line) the current as a function of voltage for the high 
dissipation limit where the harmonic oscillator goes to the 
ground state between tunneling events. 

In conclusion, we have studied fully quantized me- 
chanical motion of a single-level quantum dot coupled 
to two voltage biased electronic leads. We have shown 
that above a certain threshold voltage the energy accu- 
mulated in the mechanical subsystem increases dramat- 
ically. We have also shown that second order inelastic 
tunneling events are enough to stabilize this pumping of 
energy. Finally the current through the system was cal- 
culated and it was found that the development of the 
mechanical instability is accompanied by a dramatic in- 
crease in the current. 
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